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The CP- violating charge asymmetry in decays — ► 7r =1: 77, which is induced by the 
electromagnetic penguin operators, has been studied both in the standard model and 
in its extensions. Because of a large enhancement of the Wilson coefficients of the 
electromagnetic penguin operators in the supersymmetric extensions of the standard 
model, a significant upper bound on this charge asymmetry could be expected, and thus 
high-precision measurements of this interesting CP-violating quantity might explore 
^ ' new physics effects beyond the standard model. 
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Rare kaon decays provide a very useful laboratory both to test the standard model (SM) 
and to explore new physics beyond it [TJ El H] • It is of particular interest to study the CP- 
violating effects, which arises in such weak decays from dimension-five operators, including 
the electromagnetic and chromomagnetic penguin operators (EMO and CMO) since the CP 
violation induced by these operators is suppressed in the SM; however, it could be enhanced 
in its extensions jSJ El HI M M EE EH U2] • On the other hand, present experiments, HyperCP 
[13 , and future NA48 experiments [14| . are going to substantially improve the present limits 
on the Wilson coefficients of these operators by studying CP-violating charge asymmetries 
in charged kaon decays, such as K ± — ► (37r) =*=, K ± — > tc ± £ + £~ , as well as one-photon or two- 
photon radiative decays. It is expected that charged kaon decays could be an ideal framework 
to explore direct CP violation, or CP violation of pure AS = 1 origin [3 El 13 UH EEH EH]- 
The purpose of this paper is devoted to the analysis of the CP-violating charge asymmetry 
induced by the EMO in the two-photon radiative charged kaon decay, K — > 7r ± 77, both in 
the SM and in its possible extensions. 

The general weak effective Hamiltonian, contributed by the EMO and CMO, can be 
written as [Sj 

H cS = C+MQ+t/i) + C 7 -(/i)Q» + C^Q+t/i) + C70*)Q;0*) + H.c, (1) 
where are the Wilson coefficients and 



Qi = T7T2^L°^d R ± s R a^d L )F^, (2) 

107T 



= 77T2( s L(7u.vt a d R ± s R aa V t a d L )G^ u . (3) 

107T 

Here Qd = —1/3, and o^ v = i/2[j^, 7„]. Note that eq. (0) with complex Wilson coefficients 
could lead to new flavor structures beyond the SM, which generally depart from minimal 
flavor violation [IB]. It is easy to see that the SM structure, SU(2)l x U(1)y, will impose 
the following chiral suppression for these operators [TT] UH] : 



= \7I 



Cii-^{m d s L a^ u t a d R + m s s R <j^t a d L )G^ 



+C 12 -^(m d s L a llu d R + m s s R a fJ , l/ d L )F llL/ 



+ H.c, (4) 



with 



Cn(m w ) = ^-^^ 4(1 -x) 3 ' (5) 

_ , , x 2 (2-3x), 8x 3 + 5x 2 -7x 

Cl2{mw) = W^f x " "W {) 

where x = m^/rn^ and t a are the S'[/(3)-matrices. However, as we shall see, new flavor 
structures, for instance, from the supersymmetric extensions of the SM, allow us to avoid 
the chiral suppression for the operators in eq. (@J) [19]. 
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It has been known that the K — > 7r ± 77 transitions are dominated by long distance effects 
[201 E]- Within the framework of chiral perturbation theory [201 122], K ± — > 7r ± 77 receives 
the first non- vanishing contribution at 0(p 4 ) including both loops and anomalous and non- 
anomalous couterterms. The general 0(p 4 ) amplitude of the decay can be decomposed in 
the following way [20] 



M(K + (k)^n + (p) 1 (q l ,eMq 2 ,e 2 )) 
= e lf ,(qi)e 2u (q 2 



— VhQi -9i ■ <hF ) + 



m 



£ <?la<?2/3 



A" 



(7) 



with 



y 



k-(qi~ g 2 ) 



m 



K 



The physical region in the dimensionless variables y and z is given by 

O^M^A^M,^), 0<z<(l-r n ) 2 , 



(9) 



where r n = m^/mx and A(a, b, c) = a 2 + b 2 + c 2 — 2(ab + ac + 6c). Note that the invariant 
amplitudes A(y, z) from loops and non-anomalous counterterms, and C(y, z) from anomalous 
counterterms have to be symmetric under the interchange of qi and q 2 as required by Bose 
symmetry. In the SM, the 0(p A ) amplitude for A(y, z) has been given in Ref. (20], which is 



2ttz 



z)F 



z 



+ (1 



z)F(z) + cz 



(10) 



where oem = e 2 /47r, and \G$\ = 9.2 x 1CT 6 GeV -2 . We do not display the explicit expression 
for C(y,z) since it is irrelevant to the present discussion. F(z/r 2 ) and F(z) are generated 
from 7r and K loop diagrams respectively, which could be defined as 



F{x) 



1 4 -2 

1 — ^ arcsm 



1 + % | In 



1 - Jl-i/x 



m 



x < 4, 



a; > 4. 



(11) 



l + sjl- 4/x 

c in eq. (|10|) is from 0(p 4 ) non-anomalous local counterterms 
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[3(L 9 + L w ) + N u - N 1S - 2N : 
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(12) 



where L 9 and L 10 are couplings in the 0(p A ) strong chiral lagrangian [22] and iVj (z=14, 
15,18) are couplings in the 0(p A ) weak chiral lagrangian |22j . 

From eq. it is obvious that the 7r loop contribution, which is proportional to F(z/r 2 ) 
in eq. (jTOjl . will generate a CP invariant absorptive part. Thus if c has a non- vanishing phase, 
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the interference between these two parts will lead to the charge asymmetry in K — > 7r ± 77 
as follows [201 121 



ST 


\T(K + - 


-> 7T + 77) 


- r(A" _ - 


-> 7T 77) | 


2T 


t(k + - 


— > 7T + 77) 


+ r(A- - 


-> 7T _ 77) 


ST = 


\T(K+ - 


-> 7T + 77) 


- r(/r - 


-> 7T"77)| 



(13) 



= '""'^rf' 4 ' /T'^-- VCWJM - 1 - * (14) 

Thus information about the imaginary part of c is relevant to the estimate of ST. It is clear 
that Lg and L10 cannot contribute to it. In the SM, Imc is therefore governed by the weak 
coupling combination N u — Ni 5 — 2Ni$ due to the Cabbibo-Kobayashi-Maskawa (CKM) 
phase j2Zj- This point has been explored in the past literature, and some small values of the 
charge asymmetry in K —> 7r ± 77 have been obtained in Refs. [21] and (SlESj, respectively; 

SM 

< 10- 3 j2H, (15) 

and 

/rn\SM 

[^j < 10- 4 m (is) 

We would like to give some remarks here. First, eq. (j!5|) is the updated version of the 
result given in Ref. [201, and a vanishing imaginary part of A^ 14 — A^ 15 — 2iV 18 has been 
predicted by the authors of Ref. [28J using l/N c analysis. Second, in deriving eq. (fTHjl. 
the 0(p 6 ) unitarity corrections |23 EE] from the physical K ± — > 3tt vertex to K — > 7r ± 77 
have been taken into account. As an order-of- magnitude estimate for the charge asymmetry 
111 — > 7r ± 77 induced by the EMO, in the present paper we do not consider the 0(p 6 ) 
unitarity contributions since, as pointed out in Ref. 1213 , ^ does not alter significantly the 
value of the charge asymmetry obtained at 0(p A ) in chiral perturbation theory. 

Let us now try to delve into the analysis of the CP-violating charge asymmetry in K ± —>■ 
7r ± 77 arising from the EMO formulated in the general effective Hamiltonian of eq. It 
will be shown below that, this asymmetry is very small in the SM, however, it could be 
significant in the supersymmetric scenarios beyond the SM by comparison with those given 
in eqs. (I15|) and Using the same way given in Ref. 0, we first construct the effective 
lagrangian that represents the EMO. In fact there are many possible chiral realization of 
these operators. To the purpose of this paper, the leading order 0{p A ) realizations of the 
EMO, which are relevant to K — > Tc-f-f transitions, could be expressed as 

C = a 1 C 1 + a 2 C 2 , (17) 
£ i = ^T2-C±(\U(F Lflu + IPF^U) ± KFl»u + U^F Rpu U)U^)F^ + H.c, (18) 
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and 



eQo 
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■C±(\UL„L V ± \L^L V U ] )F^ + H.c, 



(19) 



where a% and a 2 are unknown coupling constants, L M = iWD^U, and (A)y = f^i^i- We 
use the standard notation in chiral perturbation theory [2H], Fhp,v = Fr^ = eQF^, D^U = 
d^U — ie[Q, U]Afj,, and Q = diag(2, —1, — 1)/3. U is a unitary 3x3 matrix with det U—l, 
which collects the Goldstone meson fields (n, K, and rj) as follows 

U = exp(iV2Q/f v ), 



1 V2 + VE 



$ = 7! A>a(x) 



V 



7T 



7T 



+ 



K° 



K + \ 



'V5 + >/6 



.228 
>/6 / 



(20) 



where the A a 's are 3x3 Gell-Mann matrices and f n ~ 93 MeV. Note that our result [eq. (jl7j) ] 
is not contrary to that in Ref. [2E] since it does not give any contribution to the imaginary 
part of the weak coupling combination N u — N 15 — 2iV 18 . However, it is easy to see that 
the effective lagrangian [eq. (|T7jl ] will give new contributions to decays K + — > 7r + 77 and 
K~ — > 7r~77, respectively, one can therefore expect that the new structures in eqs. (|T£|) and 
|TT)|) with the complex Wilson coefficients will induce a possible charge asymmetry in 

— » 7r ± 77 decays. Consequently, the interference between the amplitude for if* 
from eq. (|T7jl and the absorptive part in eq. (fTU|) will give 



7T =I: 77 



5T 



EMO 



|r(K+ -»■ tt+77) - r(K~ -> tt - 77)| 



EMO 



EM 



Gslm 



A' 



3 • 2 V/. 2 

-(l-r*) 



ImC 7 ^ -aj - f/ 2 



cte A 1/2 (l, 2, r 2 ,)^ - 1 - z)z lmF(z/r 2 n ). 



(21) 



The first observation of K + — ► 7r 
and the branching ratio of the decay has been measured j^U EDI 



"77 has been reported by BNL E787 Collaboration [2T] . 



Br (i^ 



7T 



"77) 



'1.1 ±0.3 ±0.1) x 10" 



(22) 



Thus we have 




EMO 



(2.4 ±0.7) x 10 e 



ImC+ [~a x -a 2 



(23) 



Our next task is to evaluate the magnitude of ImC:}~ and (2/3cti — 02) to check whether 



it is possible to get a significant charge asymmetry of K — > 7r ± 77 from eq. ([23)1 . However, 
since unknown constants a\ and a 2 in eqs. (JTHj) and ()19p are related to the low energy chiral 
dynamics, at the present we have no model-independent way to give a reliable determination 
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of them. In the following we will estimate them using naive dimensional analysis, within the 
chiral quark model, and employing lattice calculation, respectively. 

Naive dimensional analysis As the order-of-magnitude estimate, using naive dimensional 
analysis [31] . we can obtain 

a\ ~ a 2 ~ fn-^-, (24) 
with A x = 47r/ 7r as the chiral symmetry spontaneously breaking scale. So we get 

~ ai ~ a 2 ~ — , (25) 
An 

and 

EMO 

~ 1.8 x 10 4 GeV |ImC+|. (26) 





The chiral quark model The chiral quark model has been extensively used to study low 
energy hadronic physics involving strong and weak interactions [B^l ESI OH OH] ■ in order to 
study the direct CP violation in decays K — > tt£ + £~ {£ = e, /i), this model has been employed 
by the authors of Ref. [T2j to evaluate the bosonization of the EMO, which corresponds to 
the a 2 part in eq. (fTTj) . In the same way as that in Ref. one can obtain the effective 
lagrangian corresponding to the a\ part in eq. (JT7J). This leads to 

01 = ^ a2 = i^- (27) 

where the constituent quark mass Mq could be set about 0.3 GeV [T3]. Thus we get 

;i.6 ± 0.4) x 10 4 GeV |ImC+|. (28) 



01 \ 



\2T ) 

Lattice calculation So far there is no direct lattice calculation on a% and a 2 . However, 
note that the a 2 part in eq. (JT7J) will also contribute to transition K — > tij* — > tt£ + £~ , and 
the first lattice calculation of the matrix element of the EMO, (ir°\Q^\K ), has been done 
in Ref. [HE] . One therefore could determine the value of a 2 by comparing the result from eq. 
([17]) and that from the lattice calculation. In general, the matrix element of the EMO can 
be parametrized in terms of a suitable parameter jSl EH] 



(AQ+\K°) = ^^^B T V W K F» V . (29) 
On the other hand, using the effective lagrangian in eq. ([17]) . one can get 



'w-' jST 164 ^ (30) 
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Since Bt = 1.18 ±0.09 has been found in the lattice calculation [36 , from eq. (|2T?j) and (|3U|). 

we have 

a 2 = 0.010 ±0.001 GeV. (31) 

Unfortunately, now we are not able to use the similar way to extract any information on ai 
from the lattice calculation. However, by comparing the value of a 2 in eq. (J31|) with those 
from naive dimensional analysis [in eq. ()24|) ] and the chiral quark model [in eq. (|27|)]. one 
can find that they are of the same order of magnitude. Meanwhile, a\ from eqs. ()24|) and 
(j2ZJ) are also of the same order of magnitude. This leads to the same order of magnitude 
estimates for the charge asymmetry in eqs. (|26|) and (|28|). which are from naive dimensional 
analysis and the chiral quark model, respectively. Therefore, in general, 



EMO 

2T / 



1.0 x 10 4 GeV |ImC+| (32) 



could be expected except in a fine-tuning case in which there is an accidental cancellation 
between 2/3ai and a 2 (since we cannot reliably fix the relative sign of ai and a 2 in a model- 
independent way). 

In order to go further into the charge asymmetry in eq. ()32|) . now one has to compute 
the imaginary parts of the Wilson coefficients C!r, which are related to the short distance 
physics. In the SM it is easy to get ImC+ from eqs. (JTJ) and (JU) as 

|ImC+| SM = -^(m 8 + m d )\lm\ t C 12 |, (33) 

where = VtdV^ s . Due to the smallness of \m.\ t ~ 10~ 4 , this contribution from the SM to the 
charge asymmetry in K ± — > 7r ± 77 is strongly suppressed and could be negligible. Therefore 
in the following we will turn our attention to physics beyond the SM. 

Among the possible new physics scenarios, low energy supersymmetry (SUSY) [37], rep- 
resents one of the most interesting and consistent extensions of the SM. In generic supersym- 
metric models, the large number of new particles carrying flavor quantum numbers would 
naturally lead to large effects in CP violation and flavor- changing neutral current (FCNC) 
amplitudes [3~8] . Particularly, one can generate the enhancement of at one-loop, via 
intermediate squarks and gluinos, which is due both to the strong coupling constant and 
to the removal of chirality suppression present in the SM. Full expressions for the Wilson 
coefficients generated by gluino exchange at the SUSY scale can be found in Ref. [TUl . We 
are interested here only in the contributions proportional to l/m g , which are given by 

C£ susyK) = ± 0Er); 2 ] ^susyK,), (34) 

m g 

Ct susyK) = [(*? R )2i ± OSOiJ G S usy(x ot ), (35) 

m g 

where (S^ K )ij = (Mp) iUR /m? denotes the off-diagonal entries of the (down-type) squark mass 
matrix in the super-CKM basis, x gq = m~/m/~ with rrtg being the average gluino mass and 
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rriq the average squark mass. The explicit expressions of -Fsusy(^) and Gsvsy(x) are given in 
Ref. jS], but noting that they do not depend strongly on x, it is sufficient, for our purposes, 
to approximate F S u SY (x) ~ -Fsusy(I) = 2/9 and Gsusy(^) ~ Ctsusy(I) = —5/18. In any case 
it will be easy to extend the numerology once x gq is better known. Also the determination 

and ()35|) can be improved by the renormalization 
m-g = 500 GeV, m t = 174 GeV, m b = 5 GeV, and 



of the Wilson coefficients in eqs. 
group analysis [^HHH]- Then by taking 
li = m c 



1.25 GeV, we will have 



ImCj 



SUSY 



2.4 x 10" 4 GeV" 



Im[(^ R ) 21 + («$£ 



lr;i2J 



Using the experimental upper bound on the branching ratio of Kl 
KTeV Collaboration [31], the lattice calculation [36J has given 



(36) 

e~ measured by 



Im[(* 



•D 

LR. 



'21 + 



12. 



< 1.0 x 10~ 3 (95% C.L. 



(37) 



Thus from eqs. (|3~2j). (|3*B)l. and (|3*7|) . the charge asymmetry in K ± — > 7r ± 77 induced by the 
EMO in the supersymmetric extensions of the SM could be bound as 



— , < a few x 10 , 

/ SUSY 



(3? 



which is significantly larger than the charge asymmetries given in the SM (HI E3J EH] . 

In conclusion, we have studied the CP-violating charge asymmetry induced by the elec- 
tromagnetic penguin operators in K — > 7r ± 77 transitions, and supersymmetric extensions 
of the SM may enhance the Wilson coefficients of these operators, which leads to interesting 
phenomenology in this study. It is found that the constrain imposed by experiments jSHj on 
the SUSY parameter Ihi[(5lr,)21 + (^lr)* 2 ] allows a significant upper bound on the charge 
asymmetry given in eq. (|38|). Our analysis shows that the charge asymmetry in K — > 7r ± 77 
up to 10 -3 would be a signal of new physics, and thus high-precision measurements of this 
CP-violating observable might probe interesting extensions of the SM. 
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